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Abstract. In this paper we try to establish a connection between a three- 
dimensional Lotka-Volterra dynamical system and two-dimensional topo- 
logical surgery. There are many physical phenomena exhibiting two-dimensional 
topological surgery through a 'hole driUing' process. By our connection, 
such phenomena may be modelled mathematically by the above dynamical 
system. 

1. Introduction 

Topological surgery is a technique used for changing the homeomorphism 
type of a topological space. For example, all orientable surfaces may arise from 
the 2- dimensional sphere using surgery. Further, the phenomenon of surgery 
appears in many physical processes. In Section 2 we discuss in detail the 
mathematical definition of topological surgery in dimensions one, two and three 
and we give examples of physical processes exhibiting surgery in dimensions 
one and two, such as: DNA recombination, reconnection of cosmic magnetic 
lines, cell mitosis, formation of tornadoes and whirls, falaco solitons, etc. 

In Section 3 we present the three-dimensional Lotka-Volterra dynamical 
system, given and analyzed by N.Samardzija and D.Greller in [Tl] (see system 
3.1). This system generalizes the classical two-dimensional Lotka-Volterra sys- 
tem [11112] and it depends on three parameters. As it turns out, in certain para- 
metrical regions the system exhibits chaotic solutions shaping a fractal toroidal 
bead in 3-dimensional phase space. Moreover, with a continuous change of pa- 
rameters, the chaotic attractor passes from spherical to toroidal shape via a 
'drilling' process along a slow manifold, thus performing 2-dimensional topo- 
logical surgery. This connection was first observed in |i8j and consequently 
worked in [1]. In this paper we establish this connection between the paramet- 
rical behavior of the solutions and the 2-dimensional topological surgery by 
solving the system numerically in each parametrical region and demonstrat- 
ing this 'hole drilling' process. This means that system 3.1 (more precisely, a 
one-parameter family) can model 2-dimensional surgery. 

Further, there are physical phenomena resembling 2-dimensional surgery 
through a 'drilling' process, such as formation of tornadoes and whirls, falaco 
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solitons, etc. These applications are discussed in Section 4. By our connection, 
such phenomena can be modelled by system 3.1. Finally, in Section 5 we 
discuss the vectorial dynamo equation, the solutions of which we conjecture to 
model 3-dimensional surgery. This will be the subject of future work. 

2. The topological surgery 

2.1. A homeomorphism between two m-dimensional topological spaces is a 
continuous bijection, such that the inverse is also a continuous map. In Topol- 
ogy, two homeomorphic topological spaces are considered the same. For this 
reason they may be attached together and a homeomorphism between them 
can be used as 'glue'. 

An m-manifold M is a metric space which may be covered by open sets, 
each of which is homeomorphic with M"* or the half-space M™ (in case M 
has boundary). The aim of Topology is to classify all m- manifolds up to 
homeomorphism, that is, to characterize and list all non-homeomorphic m- 
manifolds, for each m. Here we restrict our interest to orientable manifolds. 

The only connected 1-manifold is the circle S^. The connected, orientable 

2- manifolds without boundary are the 2-dimensional sphere S"^ (only the sur- 
face, without the inside), the torus (without the inside), the torus with two 
holes, with three holes, etc. In three dimensions the basic 3-manifold is the 

3- dimensional sphere S^, which may be viewed as the Euclidean space with 
all points at infinity compressed into one single point. (Analogous views are 
true for spheres of all dimensions). Other examples of 3-manifolds are the lens 
spaces, homology 3-spheres, etc. 

In order to create new topological spaces out of known ones we use the 
following technique. 

Definition 1. A topological surgery on an m-manifold M is the topological 
procedure of creating a new m-manifold M' by removing the product set S*" x 
jjm-n q£ n-dimensional sphere with the (m— n)-dimensional ball, taking the 
closure of the remaining manifold with boundary and replacing S"" x ]J^~^ by 
D""*"^ X S""^"""^, using a 'gluing' homeomorphism along the common boundary 
gn X ^m-n-i^ that is: 

M' = (M \ X D'"-") U5nx5— n-i X S"^-"-! 

An excellent reference on the subject is ^U\. In this paper we shall use the 
term 'm-dimensional topological surgery' with no risk of confusion about the 
sub-dimension n. One-dimensional surgery, for example, means that we remove 
from one circle or from two circles two arcs x {S^ is two points) and we 
reconnect the four boundary points x 5° in a different way, obtaining in the 
end two or one circle depending on the type of reconnection. View Figure 1 
(cf. p. 247 (TUj). 
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to form S"*" X 



Figure 1. Surgery in dimension 1 
Source : D. Rolfsen, Knots and Links, Publish or Perish Inc. (1976), p. 247. 



In this paper we are mainly interested in 2-dimensional surgery. From Def- 
inition 1 we have: 

M' = (M\SO X D^) X 

That is, we remove two discs x from the surface M and we replace them 
by a tube x S^, which gets attached along the common boundary S'^ x S^. 
This operation changes the homeomorphism type of the initial surface. In fact, 
according to the surface classification theorem, every connected, orientable 2- 
dimensional surface without boundary arises from the sphere S"^ by repeated 
surgeries. Note that, every time the above process is performed a new 'handle' 
is added on the surface. 

For example, for M = S"^, the 2-dimensional sphere, the process is illustrated 
in Figure 2. Two points are specified first on the sphere, which are the centers 
of the discs to be removed. After the removal of the two discs x a tube 
X is attached along the common boundary S'^ x S^, resulting in a torus. 
This is a different 2-dimensional surface, arising from the sphere via the surgery 
technique. Note that, if the tube were attached on the two-punctured sphere 
externally, the result would be a homeomorphic torus. As another example, 
if the initial surface M is two copies of the sphere S"^ and the two discs to be 
removed are one on each sphere, replacing the two discs by a tube yields one 
copy of the sphere 5*^. 
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Three-dimensional topological surgery is much harder to visualize than 2- 
dimensional surgery. Letting m = 3 and n = 1 in Definition 1, it means that 
we remove an embedded solid torus (a torus with the inside) x from 
a 3-manifold M and we replace it by another solid torus x with the 
factors reversed. (To see that x is a solid torus consider a disc in 
every point of the circle S^.) The replacing of the solid torus takes place by 
attaching along the common boundary x (which is a torus) via a gluing 
homeomorphism : 



The result is a new 3-manifold M'. Taking the 3-sphere = R'^ U {oo} as our 
starting manifold, by known results in Topology, every 3-dimensional closed 
(compact, without boundary), connected, orientable manifold may arise as the 
result of surgery in along a knot or link (these are embedded circles or copies 
of circles in three-space). 

2.2. Surgery is not just a technique used in Topology for changing the home- 
omorphism type of a topological space. It also happens in nature in diverse 
processes for ensuring new results. For example, 1-dimensional topological 
surgery happens in DNA recombination for changing the genetic sequence. 
View Figure 3 (cf. Wikipedia). It also happens in the magnetic reconnection, 
the phenomenon whereby cosmic magnetic field lines from different magnetic 
domains are spliced to one another, changing their patterns of connectivity 
with respect to the sources, see Figure 4 (cf. [2]). 



Figure 3. Crossing over of chromosomes in DNA recombination 

Source : Wikipedia 



Two-dimensional topological surgery happens, for example, in the biological 
process of mitosis, where a cell splits into two new cells. View Figure 5 (cf. p. 
395 [6j). A similar picture can be found in the mechanism of gene transfer in 
bacteria (cf. p. 471 [3]). More examples are presented in Section 4. 

Finally, a possible physical application of 3-dimensional topological surgery 
is discussed in Section 4. 



(2.1) 



M' = {M\S^x D2) U51XS1 D"^ xS 
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Figure 4. The reconnection of cosmic magnetic lines 

Source : R.B. Dahlburg, S.K. Antiochos Reconnection of Antiparallel Magnetic Flux Tubes^ 
J. Geophysical Research 100, No. A9 (1995) 16991-16998. 



EARLY ANAPHASE 



LATE ANAPHASE 



TELOPHASE 



INTERPHASE 




Ccntromejes have uncoupled and 
begun moving toward opposite 
poles of spindle. 



The 2 sets of new single-stranded 
chromo!>ome.s nearing respective 
poles. 

Cytokinesis beginning. 




New nuclear membranes forming, 
Chromo,somes become longer, thinner, 

and less distinct. 
Nucleolus reappearing. 
Centrioles replicated. 
Cytokinesis nearly complete. 



Nuclear membranes complete. 
Chromosomes no longer visible. 
Cytokinesis complete. 



Figure 5. The process of mitosis is an example of 2D surgery 
Source : W.T. Keeton, C.H. McFadden, Elements of Biological Science, W.W. Norton 
Company Inc., 3rd edition (1983), p. 395. 



3. Connecting a 3-dimensional Lotka-Volterra model with 
2-dimensional topological surgery 



3.1. In [TT] Greller and Samardziya study tlie beliavior of tlie following dy- 
namical system that generalizes the classical Lotka-Volterra problem [9l [12] 
into three dimensions. We shall present the analysis of [H] using the same no- 
tations and some of their figures. The system is a two-predator and one-prey 
model, where the predators Y, Z do not interact directly with one another but 
compete for prey X. Apart from a population model, the system may also 
serve as a biological model and a chemical model (cf. [TT]). 



f ^ = X - XY + CX^ - AZX^ ) 

at 



(3.1) 



§ = -BZ + AZX'^ 



> A,B,C>0 
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The parameters A,B,C are analyzed in order to determine bifurcation prop- 
erties of the system. In particular, the phenomenon of chaos evolving on a 
fractal torus is observed and analyzed by analytically derived two-dimensional 
stability diagrams. As parameters A,B,C affect the dynamics of constituents 
X, Y, Z, the authors were able to determine conditions for which the ecosys- 
tem of the three species results in steady, periodic or chaotic behavior. More 
precisely, the authors derive the following five steady state solutions for the 
system: 



Ssl 











SsA — 1 


T 


i Ss5 — 










Since X, Y, Z are populations, the negative steady states Ss4, Ss5 are not 
taken into consideration. Figure 6 (cf. Fig. 1 in [H]) illustrates the positions 
of steady states Ssi, Ss2 and Ss2 in P^, where denotes the set of points in 
with all coordinates > 0. 




Figure 6. Repelling steady state points Ss2 and Ss3- 

Source : N. Samardzija, L. Greller, Explosive route to chaos through a fractal torus in a 
generalized Lotka-Volterra Model, Bull Math. Biology 50, No. 5 (1988), Fig. 1. 
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Let, now, J{S) be the Jacobian of system 3.1 evaluated at S" G M'^ and 
let the sets r{J(S')} and W{J{S)} to be, respectively, the eigenvalues and 
corresponding associated eigenvectors of J{S). Using the sets of eigenvalues 
and eigenvectors given below, the authors characterize the local behavior of the 
dynamical system around the points Ssi,Ss2, Sss using the Hartman-Grobman 
(or linearization) Theorem. 



r{j{Ssi)} = {i,-i,~B}- w{j{s,^)} 



1 

















1 
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r{ J(S,2)} = {A~B,{C+ v/(C-2)2-8)/2, (C - V(C - 2)2 - 8)/2} 



W{JiSs2)} 



1 

{C + 1)/{A-B) 

B+C-A+{C+l)/(B~A) 
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1 



C-a/(C-2)2-8 
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1 



C+v/(C-2)2-8 



2 





IB -1 + ^l-8B{l + C,/BjA) -1 
A~^' 2 ' 



1-8B{1 + C^/B/A) 



W{J{S,s)} 



2VAB{1+Cy 




1 





-^1-8B{1+Cy/B/A) 



2B 



1 





1-8B{1+Cy/B/A) 



2B 



Since 1 > and —1,—B < 0, Sgi is a saddle point for all A,B,C. By 
inspecting the variation of the eigenvalues for Ss2 and Ss3 in relation to A,B, C, 
the authors show that the behavior of Ss2 and Sg^ can be characterized by 
only two parameters, namely C and B/A, as shown in the stability diagrams 
in Figure 7 (cf. Figs. 4, 5 and 6 in pjj). Before examining the behavior of 
the system in each different parametrical region, the authors analyze it into 
three 2-dimensional subsystems, concluding that all solutions initiated in 
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will remain in for all times. Furthermore, they show that system 3.1) is 
always bounded in P^. 



c 





a Id g 

1 C.2(1+/2) 
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case 1/8 < 5 — stability diagram 
C 1 




la- 1 d g 

' \ 1 C = 2(l + J2) 

b la h 

^ C.{-1+1/8B) VaTB 

h 





case l/(8(3 + y2))<B< 1/8— stability diagram 
c = (-1+1 /8B) Va7b A 1 

C = 2(U,/2) / 




1 

k m n 

1 p B/A 


o 





case 05:^^1/(8(3-1-^2)) — stability diagram 



Figure 7. Stability diagrams 

Source : N. Samardzija, L. Greller, Explosive route to chaos through a fractal torus in a 
generalized Lotka-Volterra Model, BuU. Math. Biology 50, No. 5 (1988), Figs. 4, 5 and 6. 
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Parametrical regions a,b,c,d,e, f,j,k,l,m,o and r are stable steady state 
regions, while i and p are periodic/stable steady state regions. Interesting 
behavior arises in the parametrical regions g, h, n and q where Ss2 and Sg^ are 
an inward instable vortex and an outward stable vortex, respectively (view 
Figure 6). Let Ai,A2,A3 denote the eigenvalues of either 5*^2 or Ssz- In both 
cases they must satisfy the conditions Ai e M and A2, A3 e C with A3 = Ag, 
the conjugate of A2. The eigenvalues of Ss2 must further satisfy Ai < and 
Re{\2) — Re{\^) > 0, while the eigenvalues of 5*53 must further satisfy Ai > 
and Re{\2) = ReiXs) < 0. 

Note now that the point Ss2 as well as the eigenvectors corresponding to its 
two complex eigenvalues, all lie in the xy-plane. On the other hand, the point 
Ss3 and also the eigenvectors corresponding to its two complex eigenvalues 
all lie in the x2;-planc. Finally, the eigenvector corresponding to the real 
eigenvalue of Sgs breaks out of the a;2;-plane and redirects the flow towards 
Ss2- Thus, since we have two repelling steady points and since there are no 
sinks in and the system is bounded, the formation of chaotic/periodic type 
trajectories is expected in regions g, h, n and q. 

Setting B/A = 1 and equating the right side of system 3.1 to zero, one finds 
as solution the one-dimensional manifold 

L = {{X,Y,Zy, X = 1,Y + AZ = 1 + C} 

that passes through Ss2 and Sgs- Since all points on L are singular, there is no 
motion along it. If one moves in parameter space to the left of B/A — 1, the 
manifold L exerts flow toward the inward stable vortex Ss3 and if one moves 
to the right of B/A = 1, the manifold L exerts motion toward the inward 
unstable vortex Ss2- Consequently, stable solutions are generated left of and 
including the line B/A — 1, while the chaotic/periodic regions g,h,n and q 
that we are interested in appear on the right of the hue B/A — 1. 

Using CRAY FORTRAN, the authors solved system 3.1 numerically in 
chaotic/periodic regions g,h,n and q (where B/A > 1) and showed that the 
trajectories are quickly trapped into an attractor between 5*^2 and Sss- Since 
the inter-arrival times exhibit stochastic behavior, the authors state that the 
attractor is non-periodic, thus being chaotic. The line B/A = 1 exhibits the 
slowest dynamic on the attractor, so L is referred to as the 'slow manifold', 
while the 'fast manifold' is the outer shell of the attractor. As parameter 
A changes, the slow manifold L opens into a hole and the chaotic attractor 
evolves into a conventional toroidal shape. 

3.2. We arc now ready to make the connection of the two areas. Using MAT- 
LAB ODE Solver wc solved system 3.1 numerically and wc represented the 
attractor in the chaotic/periodic regions g, h,n and q. We too observed that, 
as parameter A is increased, the slow manifold L opens into a hole. Figure 8 
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illustrates instances of this process in chaotic region q for B = 0.0145, C = 5.5 
and for A = 0.01305,0.01335,0.01365,0.01395 and 0.01425. Once the trajec- 
tory is trapped into the attractor, the outer shell and the line L (slow manifold) 
are drawn. In the first instance the trajectory is creating the 2-sphere and the 
slow manifold L, which links the two steady state points Ss2 and Sss- The line 
L 'selects' the centers of the two discs x that will be removed through 
the surgery process. By increasing parameter A, the hole is drilled and our 
attractor passes from a spherical to a toroidal shape, so 2-dimensional surgery 
is performed. The same shape transition also happens in regions g, h and 
n for appropriate choices of parameters. We chose chaotic region q because 
it demonstrates our point more clearly. Therefore, system 3.1 may perform 
2-dimensional surgery in the above parametrical regions. 

The use of different numerical methods may affect the shape of the attractor 
but the 'hole drilling' process is a always a common feature. 

Conclusion: The (one-parameter family of the) three dimensional Lotka-Volterra 
system 3.1 models the process of 2-dimensional topological surgery. On the 
other hand, instead of viewing surgery as an abstract topological process, we 
may now view it as a property of an one-parameter dynamical system. 

4. Applications 

Many physical phenomena start by selecting two points (poles) and by join- 
ing them through a helicoid line. This process resembles surgery, regardless 
of the topological surface the two points belong to. Moreover, 'drilling' with 
coiling is naturally chosen in various physical processes, since it is, from the 
energy point of view, the most economical way to open a hole. Thus, as a con- 
sequence of our connection, every physical phenomenon exhibiting the above 
process can be modelled by the three-dimensional Lotka-Volterra system 3.1. 

4.1. A phenomenon exhibiting the above process is the creation of Falaco 
solitons. These are singular surfaces connected by means of a stabilizing invis- 
ible thread or string. For a full exposition see [7]. In Figure 9 (cf. pp. 147, 150 
[7j) one can observe three pairs of these topological structures in a swimming 
pool. 

The invisible thread connecting the two singular surfaces (poles) of Falaco 
solitons corresponds to our slow manifold L. The hole drilling process seen in 
Figure 9 is exactly analogous to the one observed in Section 3 (recall Figure 8). 
Starting by the two poles, and by coiling the joining line, the hole becomes 
bigger and bigger until surgery is completed. In Figure 9 one can see that each 
soliton pair contains one inward instable vortex and one outward stable vortex, 
just like our two steady points 5*^2 and Ss3 of system 3.1 (recall Figure 6). 

As mentioned in p. 152 "a Falaco soliton is a topological object that ap- 
pears at all scales, so the concept of Falaco solitons is a universal phenomenon. 
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Figure 8. Opening of the slow manifold L into a hole as pa- 
rameter A increases. 



At the microscopic level, the method offers a view of forming spin pairs 
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FALACO SOLITONS 




Figure 9. Falaco Solitons in a swimming pool 

Source : R.M. Kichn, N on- Equilibrium Systems and Irreversible Processes^ Adventures in 
Applied Topology 1, Non Equilibrium Thermodynamics, University of Houston Copyright 

CSDC. INC, (2004), pp. 147, 150. 



is different from Cooper pairs and could offer insight to superconductivity. At 
the level of cosmology, the concept of Falaco solitons can be found in solutions 
to the Navier-Stokes equations in a rotating frame of reference. Based on 
the experimental creation of Falaco solitons in a swimming pool, it has been 
conjectured that M31 and the Milky Way galaxies could be connected by a 
topological thread". View Figure 10. 




( OSMK" STRINGS FROM Ht^BBLE 



Figure 10. The galaxies M31 and Milky Way 



4.2. Another physical phenomenon which can be modelled by system 3.1 is 
the formation of tornadoes. As Figure 11 demonstrates, if certain meteoro- 
logical conditions are met, funnel-shaped clouds start descending toward the 
ground. Once they reach it, they become tornados. This procedure can be 
viewed as a two-dimensional surgery. First the poles are chosen, one on the 
tip of the cloud and the other on the ground. Then, starting from the first 
point, the wind revolves in a helicoidal motion toward the second point, until 
the hole is drilled and, thus, 2-dimensional surgery is performed. Note that, if 
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one deforms the water surface of falaco solitons to a sphere, the visual corre- 
spondence to tornado structures between the earth surface and the tropopause 
is remarkable. 




Figure 1 1 . Funnel clouds drilling and tornado formation 

4.3. Another physical phenomenon completely analogous to the formation of 
tornadoes is the formation of whirls in the sea. 

4.4. Finally, a phenomena which is very similar to the previous ones is a 
wormhole. In physics, a wormhole is a hypothetical topological feature of 
space-time that is basically a 'shortcut' through space and time. It has at least 
two mouths which are connected to a single throat or tube. View Figure 12. 
This structure was presented by Wheeler in 1955 but was considered to be 
unattainable in a practical sense. While there is no observational evidence for 
wormholes, space-times containing wormholes are known to be valid solutions 
in general relativity. 




Figure 12. A representation of a wormhole 

5. 3-DIMENSIONAL TOPOLOGICAL SURGERY AND DYNAMO EQUATIONS 

Recall now fl2.ll) that describes 3-dimensional topological surgery. Letting 
M = 5*^, in order to cut out a solid torus from we need to know the following 
property: the complement of a solid torus Vi in is another solid torus V2, 
that contains the point at infinity, and one can write S*^ = Vi V2. The Hopf 
fibration of S^, illustrated in Figure 13, gives an idea of this splitting of S^. 
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The surgery process starts and ends with two 3-dimensional manifolds than 
can be represented by two tori. The parallels of the torus of become the 
new meridians of the torus of our new 3-manifold. 




On the other hand, similar toroidal structure can be found in dynamo theory, 
see for example [3]. More precisely, when conducting fluid flows across an 
existing magnetic field, electric currents are induced which, in turn, create 
another magnetic field. When this new magnetic field reinforces the original 
magnetic field a self-sustained dynamo is created. This is called the "Dynamo 
Theory" and it explains how the magnetic field of dense objects, such as the 
white dwarves, the sun or the earth, is sustained. 

In order to simulate the Earth's magnetic field, Glatzmaier and Roberts 
solved the dynamo equation numerically [4J. Their result is precisely a three 
dimensional manifold represented by two tori, view Figure 14. 

We believe that dynamo theory is the key to our next step: modelling 3- 
dimensional surgery. This will be the subject of future work. 
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Figure 14. Numerical solution of the dynamo equation 
Source : G.A. Glatzmaicr, P.H. Roberts, Three-dimensional self-consistent computer 
simulation of the geomagnetic field reversal^ Nature 377 (1995). 
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